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SUMMARY 

A  description  of  adaptive  beamforming  in  the  presence  of  linear  constraints 
is  given.  It  is  proved  that  if  the  number  of  constraints  does  not  exceed  a 
stated  maximum  then  the  adaptive  process  determines  one  unique  beam  pattern. 

The  proof  is  first  written  out  explicitly  using  elementary  algebra,  and  second 
in  matrix  form.  The  matrix  form  has  the  advantage  of  compactness  and  makes  pos¬ 
sible  the  writing  of  a  simple  formula  for  the  weights  which  themselves  define 
the  beam  pattern. 


Copyright 

© 

Controller  HMSO  London 
1980 


This  Memorandum  expresses  the  opinions  of  the  author  and  does  not  necessarily 
represent  the  official  view  of  the  Royal  Aircraft  Establishment. 


2 


LIST  OF  CONTENTS 

Page 

1  INTRODUCTION  3 

2  BASIC  FORMULAE  3 

3  THE  PRINCIPLE  OF  ADAPTIVE  BEAMFORMING  5 

4  EXTENSION  TO  MULTIPLE  CONSTRAINTS  AND  PROOF  THAT  ABF  HAS  A 

UNIQUE  SOLUTION  6 

5  MATRIX  FORMULATION  OF  THE  PROOF  10 

6  CONCLUDING  REMARKS  AND  CHOICE  OF  CONSTRAINTS  13 

Acknowledgments  1 5 

Report  documentation  page  inside  back  cover 


f  Accession  For 

NTIS  GRAAI 

If 

DTIC  TAB 

Unannounced 

n 

Justification— 

By - - — - - 

Distribution/  _ 

Availability  Codoa_ 
Avail  end/or 
DiSt  I  Special 


A 


R-N  135 


3 


I  INTRODUCTION 

Adaptive  beamforming  (ABF)  is  a  method  of  signal  processing,  explained  in 
section  3  below.  It  is  receiving  considerable  attention  at  the  moment  with  a 
view  to  its  possible  applications  to  various  types  of  sensor  arrays  and  especially 
passive  directional  sonar.  Over  recent  years  a  quite  voluminous  literature  has 
accumulated  concentrated  upon  the  numerical  techniques  needed  to  render  the  pro¬ 
cess  practicable.  However,  none  of  the  papers  seen  by  this  writer  has  addressed 
the  question  of  the  relationship  between  the  constraints  and  the  uniqueness  of 
the  solution,  or  of  the  number  of  constraints  to  which  it  may  be  made  subject. 

It  appears  to  be  worthwhile  to  make  good  any  such  gap,  and  this  Memorandum  does 
so,  for  the  case  of  linear  constraints,  by  the  use  of  elementary  algebra  only. 

The  proof  is  written  out  explicitly  in  section  4,  followed  by  its  briefer  trans¬ 
lation  into  matrix  language  in  section  5  where  the  formal  solution  is  also  given. 
Before  doing  this  the  basic  formulae  of  the  subject  will  be  derived  and  an 
explanation  of  ABF  provided. 


BASIC  FORMULAE 


Consider  an  array  of  n  sensors,  the  complex  output  (typically  the  Fourier 
transform  of  the  time  series  derived  by  periodically  sampling  an  actual  signal) 
of  the  kth  sensor  being  Xj^  ,  k  ■  1,  2,  ...»  n.  In  the  associated  signal 
processing  the  complex  conjugates  of  the  x^  ,  denoted  by  x£  ,  are  multiplied  by 
complex  weights  w^  appropriate  to  a  particular  ’look’  direction  (or  beam),  and 
sumned  to  give  the  output 

n 


X  ^ 

k-1 


(1) 


corresponding  to  which  the  output  power  in  the  beam  which  has  been  formed  is 

n  n 


|y|2  -  yy*  -  X  xk”k  X  x‘"‘ 

k-1  f-1 


co 


SB 

OS 


.5, 


xIx*wkwt 


(2) 


What  the  processer  actually  delivers  as  an  output  for  display  is  an  average 
of  the  power  over  some  period  of  time.  Such  averages  are  statistical  expectations. 
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denoted  for  convenience  by  E  ,  so  that 


n 

e( |y I 2>  -  ^  E(xkxi)wkw* 

kfl”  1 


(3) 


where 


E<*t*, 


(4) 


p(x*,x^)  being  the  probability  density  function,  and  the  integrations  being 
taken  over  all  values  of  for  *11  directions.  Then,  taking  the  complex 

conjugate  of  (4) , 


[e(x£v]  *  JJX kx>(vx*)dVxS 


(5) 


where  p  ,  which  is  just  a  real  number  between  0  and  I,  remains  unaltered  during 
conjugation.  But  the  probability  distributions  of  and  x£  must  be 

identical  (if  xfc  is  specified  so  is  its  conjugate  x£  ,  and  conversely), 
ie  p(x£,x^)  -  p(x^,x*)  ,  which  substituted  into  (5)  gives 


which  by  (4) 


[E<xiiv]  -  //‘kv(v*;)dv> 

■  E<vt>  ■  e(xK> 


* 

i  * 


Hence,  writing  for  brevity  e^  =  E(x£x^)  ,  it  follows  that 


etk  "  ekt  ’ 


and  substitution  into  (3)  gives  for  the  expected  power 


E( !y | 2) 


L  wJ 


k,l-1 


(6) 


(7) 


R-N  135 


5 


3  THE  PRINCIPLE  OF  ADAPTIVE  BEAMFORMING 

As  usually  formulated  the  objective  of  ABF  is  that  of  minimising  the 
power  given  by  the  expression  (7)  above  while  maintaining  unit  response  in  one 
specified  direction.  This  last  requirement  may  be  stated  by  means  of  a  linear 
relationship  between  the  w^  as  can  be  seen  by  the  following  simple  example :- 
suppose  that,  for  an  incident  wave  of  given  direction  and  of  frequency  f  ,  a 
time  delay  t^  occurs  between  the  first  and  kth  sensors  for  corresponding  points 
in  the  wave  front,  then  apart  from  a  proportionality  factor  (defining  the  ampli¬ 
tude  of  the  wave  and  its  phase  at  the  first  sensor)  x^  -  e  k  .  Substitution 
of  this  in  (1)  shows  that,  for  unit  response  in  the  given  direction,  the  w^ 
must  satisfy 

-2irift, 


which  is  the  linear  relationship  referred  to.  For  conventional  beamforming  (CBF) 
w^  -  T(k)e  k  where  T(k)  ,  called  the  taper  function,  satisfies 

X"1  T(k)  ■  1  consistently  with  (8). 


The  freedom  of  choice  in  T(k)  makes  possible  the  variety  of  beam 
patterns  available  to  CBF.  There  are  serious  limitations  however.  Thus,  taking 
for  example  the  case  of  a  linear  array,  if  diffraction  secondaries  (repeats 
of  the  main  beam)  are  to  be  excluded  from  look  directions  within  ±90°  of  broad¬ 
side  the  hydrophone  spacing  nust  be  less  than  a  wavelength.  Hence,  for  a  res¬ 
tricted  number  of  hydrophones  (and  in  practice  such  restriction  is  clearly 
unavoidable),  the  array  length  will  be  correspondingly  limited.  Under  these 
circumstances  any  attempt  to  improve  sensitivity  by  narrowing  the  main  beam 
results  in  increased  side  lobes  -  and  conversely.  In  practice  there  always  will 
be  side  lobes,  and  the  further  the  main  beam  is  tilted  away  from  broadside  the 
wider  it  is,  and  the  larger  some  of  the  side  lobes  become.  If  then,  while  an 
attempt  were  being  made  to  detect  a  signal  of  interest  in  a  certain  direction 
by  steering  the  main  beam  there,  a  sufficiently  powerful  noise  was  being  picked 
up  by  such  a  side  lobe,  the  whole  operation  of  the  array  could  be  confused.  ABF 
deals  with  this  problem  by  fixing  the  response  of  the  array  in  the  direction  of 
interest,  using  a  linear  constraint  on  the  weighting  factors  w^  (as  illustrated 


by  (8)  above),  and  then  exploiting  the  remaining  available  variability  of  the  w^ 
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in  re-shaping  the  beam  pattern  so  as  to  minimize  the  total  energy  received  from 
all  directions. 

4  EXTENSION  TO  MULTIPLE  CONSTRAINTS  AND  PROOF  THAT  ABF  HAS  A  UNIQUE  SOLUTION 

A  generalization  of  the  above  ABF  method  to  cover  multiple  constraints  will 
now  be  considered.  It  will  be  proved  that  any  number  m  <  n  -  1  (where  n  is 
the  number  of  independent  sensors)  of  arbitrary  independent  linear  constraints 
may  be  imposed  on  the  weights  w^  ,  k  »  1,  2,  ....  n,  subject  to  which  there 
exists  a  unique  set  of  values  of  the  w^  (given  explicitly  by  formula  (28)  in 
section  5)  for  which  the  power  given  by  expression  (7)  attains  a  minimum. 

Let  there  be  m  postulated  constraints  which  define  the  array  response 
for  m  specified  directions  such  that 


n 

£ 


* 

c  .w. 
rj  J 


(9) 


being  the  prescribed  response  in  the  rth  direction. 


where  r  ■  I ,  2,  . . . ,  m; 
the  direction  being  determined  by  the  coefficient  crj  .  So  far  as  the  present 
argument  is  concerned  the  c.j  and  qr  are  quite  arbitrary,  and  do  not  need  to 
have  the  merely  typical  physical  interpretations  just  given  (vide  section  6). 
Clearly  m  <  n  -  I  since  there  are  only  n  available  w^s  in  (9),  and  if 
m  ■  n  they  would  be  completely  determined  by  equations  (9)  alone  leaving  no 
variability  for  attaining  a  minimum  of  the  expression  (7).  In  practice  m 
would  probably  be  much  smaller  than  n  to  provide  as  much  freedom  as  possible 
for  variation  of  the  Wj  in  making  the  minimum  of  E(|y|  )  as  small  as  possible. 

The  complex  conjugates  of  equations  (9),  namely 


£ 

j-1 


* 

c  .w. 
rj  J 


1,2,. ...m  <n  -  1)  (10) 


may  also  be  written. 

The  m  constraint  equations  (9)  leave  n  -  m  of  the  w.  free  to  be  varied 

J  J 

in  search  of  minima  of  E(|y|  )  .  This  observation  suggests  that  a  way  to  find 
a  minimum  would  be  to  use  the  equations  (9)  to  eliminate  m  of  the  w.  from 
the  formula  (7)  leaving  E ( | y |  )  expressed  in  terms  of  the  remaining  n  -  m 
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variables.  If  each  of  these  remaining  w.  s  were  then  arbitrarily  varied  and 

•  .  .  .  i  1 2 

the  consequential  first-order  variations  in  E ( | y j  )  set  to  zero  there  would 

result  n  -  m  linear  equations  to  determine  the  values  of  the  n  -  m  w^  s. 

These  determined  values  when  substituted  into  equations  (9)  would  provide  ra 

equations  to  give  the  values  of  the  m  w.  s  which  were  originally  eliminated. 

This  discussion  incidentally  makes  an  argument  for  the  existence  of  a  unique  set 

of  w.  giving  a  stationary  value  (not  necessarily  a  maximum  or  a  minimum,  the 
J  o 

existence  of  which  calls  for  further  enquiry)  of  E(|y|").  When  translated  into 

the  actual  algebraic  steps  implied  however  the  formulation  becomes  unattractively 

complicated.  Also,  in  practice,  no  guidance  could  be  given  as  to  which  of  the 

Wj  should  first  be  eliminated,  and  if  an  unfortunate  choice  were  made  then  badly 

conditioned  equations  might  appear  for  the  remaining  w^  leading  to  numerical 

problems. 


These  difficulties  may  be  circumvented  and  complete  symmetry  maintained 

throughout  by  treating  all  the  w.  on  an  equal  footing  from  the  start.  Then 

J  2 

setting  to  zero  the  first-order  changes  in  E(|y|  )  consequent  upon  arbitrary 
variations  in  all  the  w^  would  give  n  equations,  while  the  constraints  (9) 
give  m  equations,  thus  making  in  all  m  +  n  equations  between  the  n  variables 
Wj  .  To  preserve  the  required  freedom  of  action,  m  extra  variables,  written 
A^(r  =  l,2,...,m)  are  introduced  and  the  following  n  equations  written 


(ID 


for  j  -  I ,  2, 
n  equations 


n;  with  which,  by  taking  the  complex  conjugates,  the 


ft  ft 

e .  .w. 
ij  i 


m 

I 


*  * 
c  .  X 
rJ  r 


0 


(12) 


for  j  *  I,  2,  ...,  n  are  associated. 

The  m  +  n  equations  (9)  and  (II)  together  determine  a  unique  set  of 
values  for  the  w^(i  ■  l,2,...,n)  and  the  *r(r  "  !»2,...,m):  the  actual 
formulae  for  these  unknowns,  (27)  and  (28),  are  derived  in  section  5  where  it 
is  seen  that  the  values  of  w^  do  not  depend  in  any  way  upon  the  which  dis¬ 

appear  from  the  argument  as  soon  as  they  have  served  their  purpose  of  maintaining 
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its  symmetry.  Henceforth  in  this  section  the  w.  will  be  assumed  to  stand  for 

i  2 

this  particular  set  of  values  which  it  will  be  proved  give  a  minimum  for  E(|y|  ). 

To  this  end  vary  these  w.  by  arbitrary  amounts  h.  ,  subject  only  to  the 

constraints  (9),  to  become  w.  +  h.  (i  ■  l,2,...,n).  Let  E(|y|  )  stand  for 

the  value  of  E(|y|  )  (defined  by  equation  (7>)  when  expressed  in  terms  of  the 
2 

w^  ,  and  E(|y|  ^ts  value  when  the 

the  variation  in  E(|y|2)  is,  by  use  of 


w^  are  replaced  by  w^  +  Ik  .  Then 
(7) 


SE(|y|2) 


61 1  (Wlh|  +  "l11^  +  el2(wlh2  +  +  '*•  +  eln(wlhn  +  WnV 

+  e2](w2h*  +  w*h2)  +  e22(w2h*  +  w*h2)  +  ...  +  e^w^  +  w%2> 


+  e  , (w  h, 
nl  n  1 


e„2(w„h2 


+  w*h 
2  r 


)  + 


+  e  (w  h  + 
nn  n  n 


w*h  ) 
n  n 


e. .h.h* 
iJ  i  J 


03) 


2 

The  right  hand  side  of  (13)  expresses  SE(|y|  )  as  a  sum  of  the  first 
degree  terms  in  h^  and  h*  written  in  extenso,  and  of  the  second  degree  terms 
h.h*  under  the  summation  sign.  Concentrating  to  begin  with  upon  the  first 
degree  terms,  and  collecting  together  separately  all  those  in  h^(i  *  l,2,...,n) 
and  those  in  h*(i  *  l,2,...,n)  they  amount  to 


which  by  use  of  (6)  is  equal  to 


! 
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£4 -CH 


i-1 

n 


hj  +  . • .  + 


Ih^  +  ...  + I 


f  ..  ,1* 

/ 

ye*A  . 

jn  J  /  n 

/ 


By  substitution  from  (11)  and  (12)  the  last  expression  can  be  written 


\  /  r-»  \  \ 

,S  Cr,T  IS  Cr2>tr . (Sc"lr 

(£4 -(£4 . (£4 


which  upon  re-arrangement  in  terms  of  X^  ,  X^  becomes 


n 

n 

fiat 

n 

I 

c. .h*  +  X 
ij  J  2 

l 

C2jhj 

+  x  y 

,  * 
c  .h. 
mJ  J 

j-1 

j  =  l 

n 

n 

n 

*;£ 

cl.h.  +  X* 
ij  J  2 

l 

c*.„.  ♦  ... 

♦  x*  y 

c*.h. 
mj  J 

j-i 

j-i 

and  the 

w.  +  h. 

1  1 

are 

required  to 

satisfy 

the  constraints 

J 

Substitution  then  of  each  of  these  two  sets  of  variables  into  each  of  the  m 
equations  (9)  followed  by  subtraction  leaves  the  m  results 


c*.h. 
rJ  J 


0 


for  r  -  I ,  2, 
accordingly 


m  . 


Similarly  the 


w* 


J 


and  w*  +  h* 
J  J 


satisfy  (10),  and 


(14) 
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I 

i  =  l 


c  .h . 
H  3 


0  , 


also  for  r  =  1,  2 . m  .  Thus,  returning  to  the  expression  (14)  it  is  seen 

•  a 

that  the  coefficients  of  all  the  A  and  A  are  zero,  and  hence  the  expression 

r  r  2 

vanishes  identically  -  there  is  therefore  no  first  order  variation  in  6E(|y|  )  , 
2 

and  E ( | y |  )  has  a  stationary  value  at  the  unique  set  of  values  of 
w.(i  =  l,2,...,n)  determined  by  (9)  and  (11). 


The  only  remaining  contribution  to  the  variation  in  E(|y|  )  caused  by 


the  h.  and  hv  is  the  last  term  written  on  the  right  hand  side  of  (13).  This 
ii  % 

is  an  exact  replica  of  the  expression  (7)  for  E ( | y  j  )  with  the  w  ,w  replaced 

2iC  X 

oy  n.,  n.  .  But  |y|  by  its  very  contrivance  is  inherently  and  inevitably 
1  3  2 

positive,  and  so  also  therefore  must  be  its  expectation  E ( j y |  )  .  Thus  the 

second  order  variation  in  the  power  is  positive  whatever  the  values  of  h^  (and 

their  accompanying  h*)  so  long  as  the  weights  satisfy  the  equations  (11) 

and  both  the  w.  and  their  varied  values  w.  +  h.  satisfy  the  constraints  (9). 

In  other  words  a  unique  minimum  exists  for  the  expected  power  for  the  chosen 

look  direction,  and  this  corresponds  to  a  unique  set  of  weights  w^  with  its 

associated  beam  pattern. 


5  MATRIX  FORMULATION  OF  THE  PROOF 

This  section  is  devoted  to  a  rendering  of  the  rather  lengthy  exercise  of 
section  4  into  the  much  shorter  matrix  version.  It  incidentally  makes  possible 
the  writing  of  a  deceptively  compact  formula  (28),  for  the  optimum  values  of  the 
weights  sought  by  the  ABF  process. 

As  before  the  basic  formulae  derived  in  section  2  provide  the  starting 
point.  The  n  complex  outputs  x^(k  =  l,2,...,n)  are  written  as  a  column 
vector  x  ,  and  the  n  complex  weights  w^  similarly  as  w  .  The  letter  T 
standing  above  and  to  the  right  of  a  matrix  or  vector  denotes  its  complex 
transpose,  ie  its  transpose  with  its  elements  replaced  by  their  complex  conjugates, 

also  known  as  the  Hermitian  transpose.  Thus  the  row  vector  of  elements  x*  is 

,  X  T  k 

written  x  ,  so  that  equation  (1)  appears  at  y  ■  x  w  ,  and  the  expression  (2) 

for  the  power  as 


2 


T  T 

W  XX  w 


(15) 
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T  .  T 

The  product  xx  of  the  column  vector  x  with  the  row  vector  x  appear¬ 
ing  in  (15)  is  a  square  n  *  n  matrix  of  elements  x^x*  •  Hence  formula  (3) 
for  the  expectation  of  the  power  may  be  written 

E  (| y| 2)  =  wTE(xx^)w 

=  w^Rw  ,  (16) 

where  R  is  a  square  matrix  of  n  rows  and  n  columns  of  elements  E(x^x*)  = e 
according  to  the  definition  of  e^  given  toward  the  end  of  section  2.  R  is 
called  the  covariance  matrix  of  the  sensor  outputs  x^  . 

From  (6)  it  can  be  seen  at  once  that 

RT  =  R  ,  (17) 

which  is  the  defining  property  of  a  Hermitian  matrix,  and  is  essential  to  the 
argument . 

The  set  of  m  constraints  (9)  takes  the  form 

CTw  =  a  (18) 

where  C  =  [c j r ]  ( j  =  1 , 2 , . . . , n  :  r  =  1 ,2 . m)  is  a  matrix  of  n  rows  and  m 

columns,  and  q  =  (qr)  is  a  column  vector  of  the  m  stipulated  responses. 

The  matrix  equation 

Rw  -  CA  =  0  ,  (19) 

involving  the  extra  column  vector  A  =  (A^)  of  unknowns  A^_(r  =  |s2 . m)  as 

well  as  that  of  the  w^  ,  is  then  written.  Equations  (18)  and  (19)  are  a  pair 

of  simultaneous  equations  in  the  vectors  w  and  A  which  determine  both  of 

them,  as  will  be  demonstrated  shortly.  Before  doing  this,  however,  it  will  be 

2 

shown  that  the  w  thus  determined,  the  optimum  w^  say,  makes  the  E(|y|  )  in 
(16)  a  minimum. 

To  this  end  let  w^  be  altered  by  a  column  vector  h  =  (h^)(k  =  l,2,...,n) 

which  is  arbitrary  apart  from  w  +  h  needing  to  satisfy  the  constraint  (18). 

— 'J  2 

Then  the  consequential  change  in  E(|y|  ),  found  by  substituting  w^  +  h  and 
Wq  successively  into  the  right  hand  side  of  (16)  and  subtracting,  is 
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6E ( | y | 2 )  =  <Wq  +  h)TR(wQ  +  h)  -  WqRwq 

=  hTRwQ  +  WgRh  +  hTRh  .  (20) 

But  since  satisfies  (19)  then 

Rwq  =  CX  (21) 

with  its  Hermitian  transpose 


xV 


(22) 


T 

remembering  R  =  R  from  (17). 

Substituting  (21)  into  the  first  term  and  (22)  into  the  second  term  on  the  right 
of  (20)  gives 


i  1 2  T  T  T  T 

6E  ( |  y  |  )  =  h  CX  +  X  C  h  +  h  Rh  . 

Again,  both  wQ  +  h  and  satisfy  the  constraint  (18),  ie 

CT(w0  +  h)  =  CTw0  +  CTh  =  a 

and 

=  a  ’ 

whence  by  subtraction 

T 

C  h  =•  0 

of  which  the  Hermitian  transpose 

T 

h  C  =  0 


(23) 


(24) 


(25) 


(26) 


must  also  of  course  be  true. 

Substitution  from  (25)  and  (26)  into  respectively  the  second  and  first  terms 
on  the  right  of  (23)  then  leaves 

6E( | y ] 2 )  -  hTRh 


I 
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only,  which  is  just  the  expression  on  the  right  of  (16)  for 

2 

replaced  by  h  .  But  E ( | y j  )  is  positive  whatever  values 
2 

take.  Thus  SE(|y|“)  is  positive  and  so,  for  w  equal  to 

2 

by  (21)  and  (24),  E | y }  )  is  a  minimum. 


i  1 2 

E ( | y |  )  with  w 
its  variables  may 
the  Wq  determined 


For  completeness 
Substitution  for  w^ 


the 

from 


actual  solutions  for  and  X 

(21)  into  (24),  or  from  (19)  into 


will  now  be  derived. 
(18)  gives 


whence 


T  -1 
C  R  CA 


q  . 


X  = 


T  -1 
C  R  C 


and  this  substituted  into  (21)  gives 


R~ 1  Cq 

'  cV'c 


(27) 


(28) 


The  formula  last  written  is  the  explicit  solution  of  the  array  processing 

problem  using  ABF  .  It  is  seen,  however,  to  require  the  evaluation  of  R  '  , 

the  inverse  of  the  covariance  matrix  whose  order  equals  the  number  n  of  sensors 

in  the  array.  Such  inversions  for  large  n  are  a  notoriously  difficult  problem. 

Because  of  this,  much  of  the  effort  in  support  of  ABF  has  been  aimed  toward 

2 

methods  of  seeking  the  minimum  of  E(|y|  )  otherwise  by  finding  the  most  rapid 

2 

descent  path  across  the  surface  of  the  function  E ( [ y |  )  in  w  space.  This 
paper  demonstrates  that  for  a  given  set  of  covariances  and  defined  constraints, 
such  paths  all  lead  to  a  unique  minimum. 

6  CONCLUDING  REMARKS  AND  CHOICE  OF  CONSTRAINTS 

This  Memorandum  gives  a  rigorous  proof  of  the  existence  of  a  unique  solu¬ 
tion  to  the  minimisation  problem  of  ABF  in  the  particular  case  of  linear  con¬ 
straints  but,  also  for  this  particular  case,  generalized  as  far  as  is  possible 
to  cover  the  maximum  number  of  such  constraints.  It  gives  a  formula,  (28),  for 
the  solution  in  terms  of  the  noise  field  as  expressed  by  the  covariance  matrix  R 
and  in  terms  of  the  constraints  through  their  directions  represented  by  the 
matrix  C  and  the  vector  of  their  magnitudes  q  . 

As  pointed  out  in  the  remarks  following  equation  (9)  in  section  4,  these 
represent  but  one  possible  choice  of  meanings  for  C  and  q  ,  and  the  mathemati¬ 
cal  argument  is  unaffected  whatever  physical  interpretation  they  are  given.  Some 
examples  of  other  meanings  are  discussed  below. 
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The  practical  significance  of  the  theorem  is  that  whatever  method  be  used 
to  seek  the  minimum  it  can  lead  only  to  one  unique  solution.  Thus,  for  example, 
if  a  procedure  of  successive  approximations  were  adopted  there  would  be  no  danger 
of  convergence  toward  some  spurious  local  minimum,  as  there  exists  none  but  just 
the  one  true  minimum  to  be  sought. 

On  the  other  hand,  the  analysis  of  this  paper  has  nothing  to  say  about  the 
beamshape  produced  by  the  particular  ABF  process  studied.  It  is  consideration 
of  this,  however,  which  can  strongly  influence  the  choice  of  constraints  in  the 
first  place.  After  all,  the  whole  purpose  of  ABF  is  to  produce  some  sort  of 
'best  possible'  beam  pattern.  If  for  example  the  signal  of  interest  should  lie 
close  to,  but  not  exactly  along,  the  direction  demanded  by  a  constraint,  then 
the  ABF  process  may  inherently  attempt  to  reject  the  genuine  signal  while  con¬ 
centrating  upon  the  constrained  direction,  forming  a  very  narrow  beam  in  the 
process.  One  way  to  avoid  this  hazard  would  be  by  requiring  the  average  response 
or  responses  over  an  assigned  sector  or  sectors  to  take  prescribed  values,  and 
this,  leading  as  it  does  to  linear  constraints,  is  covered  by  the  present  work. 

The  application  of  this  approach  is  the  subject  of  a  study  by  G.J.  Lawson  of 
Radio  and  Navigation  Department,  RAE  to  be  published  in  due  course. 

A  quite  different  approach  depends  upon  the  fact  that  very  narrow  beams 
are  generally  associated  with  the  occurrence  of  very  large  positive  and  negative 
values  of  the  weights.  The  ABF  process  could  then  be  prevented  from  generating 
such  beams  if  a  constraint  of  the  form  limiting  the  sum  of  the  squares  of  the 
weights  to  a  reasonable  value  were  imposed.  Such  quadratic  constraints  have  been 
studied  by  J.  Hudson  of  Loughborough  University  of  Technology  -  they  fall  outside 
the  scope  of  this  Memorandum. 

If  yet  another  method  of  avoiding  over-narrow  beams  were  used,  namely 
requiring  that  the  rate  of  change  of  response  with  look  direction  be  limited,  then 
although  the  constraints  are  still  linear  they  then  take  the  form  of  inequalities. 
Again  the  proof  given  above  does  not  apply.  The  same  would  be  true  if,  in  order 
to  maintain  a  prescribed  beam  width,  constraints  were  imposed  which  require  the 
responses  to  exceed  a  defined  value  in  certain  directions  close  to  the  look 
direction. 

In  short,  this  paper  establishes  the  basic  uniqueness  theorem  in  the  most 
general  form  for  ABF  with  linear  constraints,  but  makes  no  contribution  to 
the  theory  when  the  constraints  are  quadratic  or  are  expressed  as  inequalities. 
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